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Abstract 


A simple approximation is developed for the total reflectance and 
transmittance of non-absorbing turbid plane-parallel layers. It is based 
on the energy conservation law and asymptotic equations for optically 
thick turbid slabs. The accuracy of the approximation is better than 
1.5% at optical thickness larger than unity. This is comparable with the 
accuracy of many numerical schemes for the solution of the radiative 


transfer equation. 
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1. Introduction 


Radiative transfer in turbid media plays an important role in 
various branches of modern optics, including bio-optics and 
liquid crystals optics, as well as oceanic and atmospheric 
research [1]. Diffuse transmittance t and reflectance r (or 
spherical albedo) are routinely measured to extract information 
on the microstructure of the medium from light scattering 
experiments [2, 3]. Thus, it is of importance to have simple 
and yet accurate equations, which can be used for their 
calculation. This allows us to avoid numerical calculations, 
involving the solution of the integro-differential radiative 
transfer equation [1]. 


It has been found that r and ¢ depend mostly on the 
optical thickness t and the asymmetry parameter g of the 
turbid medium (both for large [1] and small [4] optical 
thickness t). We exploit these findings together with the 
energy conservation law to derive simple approximations for 


the diffuse transmittance and reflectance of turbid slabs. 


The study is limited to the case of homogeneous non- 
absorbing plane-parallel turbid media. 


2. Approximations 


The diffuse reflectance r and diffuse transmittance t are defined 
by the following equations [1]: 


m/2 
r= f do sin (200) (Yo), 
0 
1 2n m/2 (1) 
Jo) = =| ay | di sin(20)R (vo, ¥, Q), 
T Jo 0 


n/2 
t= l dv sin(2¥0)T (Vo), 
0 
1 2n m/2 (2) 
To= [do [da sin20yT O 0.9) 
27 0 0 
where 19 is the light incidence angle, 1 is the observation angle, 


9 = ~2—$, Yı is the azimuth angle of the incident light beam, 
g2 is the azimuth angle of the observation direction, 


I* (o, 0, 
R(®o, 9, p) = A p) 
I* (ðo) 
I (30, 3, p) e 
T (99, 9,9) = au 
I* (ðo) 
and [*(3)) = Foto, Ho = coso. Here x Fo is the 


hemispherical incident flux. Z} (J*) is the intensity of the 
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Figure 1. The dependence of the diffuse reflectance r, the diffuse 
transmittance ¢ and the total transmittance t, on t(1 — g). The 
results shown have been obtained by numerical solution of the 
integro-differential radiative transfer equation. Curves with the 
smaller maximum values of t correspond to the smaller values of g 
(see table 1). Curves r(t*) and t,(t*) almost coincide for all cases 
given in table 1 (except at g = 0.1268, solid curve). 


transmitted (reflected) diffused light, which can be obtained 
from the solution of the integro-differential radiative transfer 
equation. 

The results of calculations of values r, t, and additionally 
the total transmittance tj = t + tg;,, Where the diffused 
transmittance of the direct light tg;, is given by 


1 
tdir = af duo Ho exp(—T/uo), (4) 
0 


are presented in figure 1. They were obtained from the 
numerical solution of the radiative transfer equation [1] in the 
framework of the discrete ordinate method [5]. 

Calculations were performed for a plane-parallel layer of 
nonabsorbing water droplets in air. It is assumed that the 
droplet size distribution f(a) (a is the radius) is given by the 
following equation: 


f(a) = Aa’ exp(—1.5a/aep), (5) 


where A is the normalization constant ( I N f(a)da = 1). The 
effective radius aep is defined as 


= ii a’ f (a) da 
ef = i a? f (a) da’ (6) 


The distribution (5) at the effective radius aef = 6 um was 
introduced by Deirmendjian [6] for the characterization of 
droplet distributions in terrestrial clouds. 

The value of aef was not constant in our calculations. It 
was varied in the range 0.06—16.0 um. Effective radii aef used 
in calculations are presented in table 1 together with asymmetry 
parameters: 


1 
g= } [ PO sinocosoao, (7) 
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Table 1. The values of a, and g used in calculations. 


dey (um) g 
0.06 0.1268 
0.23 0.7187 
0.6 0.8289 
6.0 0.8499 

16.0 0.8691 


n On Bn 

1 0.059 7.19 
2 0.356 2.24 
3 0.584 1.17 


which were calculated using the Mie theory [6] for spherical 
polydispersions of water droplets. The wavelength was equal 
to 650 nm. Note that 0 in equation (7) is the scattering angle 
and p(@) is the phase function, which gives a conditional 
probability for photons to be scattered in a direction specified 
by the angle 0. 

The value fy;, gives the diffused transmittance of the direct 
light and t, is the total transmittance. Thus, due to the energy 
conservation law for non-absorbing layers, it follows that 


r+tt+tg, = 1. (8) 


We see that data presented in figure 1 satisfy this equation. We 
also note that the value of r is the least sensitive to the values 
of ae (and thus also g), if plotted against the transport optical 
thickness t* = t(1 — g). The value of t is most sensitive to 
the value of g, if plotted against t*. 

It is known from the general theory [1] that the value of t 
depends almost exclusively on t* as t* becomes larger than a 
fixed number tř, which is a function of g. This also follows 
from figure 1. We see that this asymptotic limit is reached for 
smaller values of tř, if g is larger. In particular, it follows from 
figure | that tř is in the range 0.5-1.5 at g ~ 0.72-0.87 and it 
is around 5 at g = 0.13. On the other hand, t, is only weakly 
sensitive to g as plotted against t*. 

The most convenient function for the approximation is the 
diffuse transmittance r (see figure 1). If it is known, then we 
have from equation (8) 


t=1—-—r-—tgiy. (9) 


The value ty; can be computed from equation (4) or by using 
the approximation 


(10) 


3 
tdir = ) ane Pat, 
n=1 


where œn, Ên atn = 1, 2,3 are given in table 2. The accuracy 
of this equation as compared to integral (4) is better than 1% 
att < 4. For larger t, we have tair < 0.0055 and, therefore, 
the direct light almost does not contribute to ts. 

It is known [1, 7] that r ~ ra, where 


o —1+0.75t* 
= , 


11 
o +0.75t* am} 
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Figure 2. The errors of equation (17) with account of 
equations (16), (10) and (11) for values of ż, t, and r at aef = 6 um. 


as t*/t;* > 1. The parameter 


1 

c= 3f du u’ Kolu) (12) 
is almost insensitive to the medium microstructure and is equal 
approximately [1, 7] to 1.07. Here u = cos Ŷ and Ko(j1) is the 
escape function, which is defined via the reflection function of 
the semi-infinite nonabsorbing turbid media Rẹœ(u, uo, p) by 
the following equation [7]: 


1 
Kolu) = o.5| +2 f duo Hp Roo (HL, wo], (13) 


where 


1 27 
Ræ(H, Ho) = zÍ Ro(L, Ho, 9) dọ. (14) 
0 


We see that r,(t = 0) = 1 — o7! > 0. It follows both 
from numerical calculations and the physics of the problem, 
however, that r(t = 0) = 0 (see figure 1). Thus, the difference 

S=fqa-r (15) 
is a positive number. This is confirmed by numerical 
calculations as well. Our task is to approximate the function 
s(t*). Then the value of r can be found from equation (15) as 
r = Ta — S (see also using equation (11)). 

Clearly, we have the following limits for the function 
s(t*): s(t*) —> Oas t* — o& (because the value of r, is the 
asymptotic result for r at large thickness, see equation (11)) 
and s(t*) > 1—o7! as t* — 0 (because then the value of 
r — 0, and thus s —> ra, for small optical thickness). These 
conditions are satisfied, e.g., by the expression 

s(t*) = Aexp(—Brt*) + Cexp(—Drt*), (16) 
which is similar to equation (10), but with other constants 
involved. In particular, we have C = 1—o~—!—A. This insures 
that s(0) = 1 —o7! and, correspondingly, r(0), calculated 
from equations (15), (16) and (11), is also zero. Constants 


optical thickness 


Figure 3. The errors of equations (17) with account of 

equations (16), (10) and (11) for values of t, ts and r at various asf, 
given in table 1. Triangles with central points correspond to the 
smallest particles (0.06 um) and filled squares correspond to the 
largest particles (16 zm). Results for other sizes of particles are 
distributed between these two limiting cases. 


A = 0.047, B = 5.26 and D = 25.0 were found from the 
fitting procedure for the case a, = 6 um (see table 1). 
Final approximate expressions have the following forms: 
r=rř; = 83 t-=1—r,; t=1-—r-—tg,, (17) 
where rq, S, tdir are given by equations (11), (16) and (10) 
respectively. We also can define t4 = 1 — ra. Then it follows 
(see equation (11)) that 


1 


ig a3" 
O+ GT 


(18) 


This is the well known result [1, 7] for the asymptotic diffused 
transmittance, which is valid as t* —> oo. 

Errors of equation (17) ata. = 6 um are given in figure 2. 
They are smaller than 1.2% for arbitrary t. This is due to the 
fact that values of constants A, B, C, D were found at this 
value of aef. It follows from figure 3 that errors somewhat 
increase for other values of aef, given in table 1. In particular, 
they reach 15% for the value of r at aef = 0.06 um (triangles 
with central points in figure 3) and 13% for the value of tẹ at 
aef = 16 um (filled squares in figure 3) at small t, where, 
however, other approximations are available [4]. The errors 
are smaller than 1.5% at t > 1 for all cases studied. 


3. Conclusion 


We propose here a simple approximation (equations (17) 
with account of equations (16), (10) and (11)) for the 
calculation of the diffuse reflectance, diffuse transmittance 
and total transmittance of turbid nonabsorbing plane-parallel 
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homogeneous slabs. The accuracy of these formulae is better 
than 1.5% at t > 1 for all cases studied. Thus, the accuracy 
is better than the accuracy of measurements involved. The 
errors are also on the order of those of most numerical radiative 
transfer equation solution results. Note that the asymptotic 
theory (see, e.g., equation (11)) can be applied with the same 
accuracy only at optical thickness larger than approximately 10 
(for water clouds [8]). Our modification can be used for optical 
thickness at least an order of magnitude smaller than that. 

Formulae presented can be used to test new numerical 
schemes of the integro-differential radiative transfer equation 
solution. They are also of help in the solution of the 
inverse problem (e.g., the derivation of the transport extinction 
coefficients t*/Z for slabs with geometrical thickness L). 

All calculations were performed for a specific particle size 
distribution, given by equation (5). However, it is known [9] 
that radiative transfer characteristics do not depend very much 
on details of particle size distributions f(a), having the same 
effective radius. Therefore, our results can be also applied to 
other types of size distribution. 

Moreover, itis known [1] that equation (1 1) is valid for any 
refractive index of the particles. Therefore, we expect that our 
results also can be applied to light scattering media different 
from clouds and fogs. 
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